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We show that the newly measured branching ratios of vector charmonia (J/if), ip' and 0(3770)) 
into jP, where P stands for light pseudoscalar mesons 7r°, n, and n' , can be well understood in 
the framework of vector meson dominance (VMD) in association with the n c -n(rj') mixings due to 
the axial gluonic anomaly. These two mechanisms behave differently in J/0 and 0' —¥ jP. A 
coherent understanding of the branching ratio patterns observed in J/0(0') — > 7P can be achieved 
by self-consistently including those transition mechanisms at hadronic level. The branching ratios 
for 0(3770) — > 7P are predicted to be rather small. 
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I. INTRODUCTION 

The recent measurements of the vector charmonium radiative decays to light pseudoscalars, i.e. J/jP, ip' and 
-0(3770) — » 77T°, 777 and 777', have brought surprises and interests to us. Earlier, the CLEO Collaboration [l| renewed 
the branchingratios for J /ip — > 77r°, 777, 777', and ip' — > 777', which arc consistent with the averages from 2008 Particle 
Data Group [2J. The branching ratio upper limits for ip' — > -fir and 777 were set, which were more than one order 
of magnitude smaller than that for ip' 77/. Meanwhile, the upper limits for -0(3770) — > 7P, where P stands for 
pseudoscalar 7r°, 77 and 77', were set to be about 10~ 5 ~ 10~ 6 . The ip' radiative decays are also investigated by the 
BESIII Collaboration with the newly collected 106 million ip' events, and the results turn out to be tantalizing. It 
shows that the branching ratios for ip' — > 771"° and 777 are only at an order of 10~ 6 , which are nearly two orders of 
magnitude smaller than ip' — > 777' 0]. 

The mysterious aspects somehow arc correlated with the J / ip and ip' data. It is found that the branching ratio for 
J /ip — > 77T° is much smaller than those for J /ip — > 777 and 777' 0, E[. This could be a consequence of suppressions 
of gluon couplings to isovector currents. As a comparison, the observation in ip' — > 7P is indeed puzzling. The 
immediate question is, what drives the difference of decay patterns between J /ip and ip' . 

In the literature, the radiative decays of the vector charmonia attracted a lot of theoretical efforts. An early study 
by the QCD sum rules [|| suggested the dominance of short-distance cc annihilations. The gluon and qq transition 
matrix elements were computed by coupling the gluon fields to the pseudoscalar states with which the branching 
ratio fraction BR(J/tp — s- jif)/BR(J/ip — s- 777) was satisfactorily described. In Ref. @, the r] c mixings with the light 
pseudoscalars rj and 77' were extracted through the axial gluonic anomaly on the basis of chiral and large N c approach. 
By assuming that the partial widths of J/ip — > 777 and 777' were saturated by the f] c -T)(r)') mixing, the branching ratios 
for J/iJj — > "frj and 777' were accounted for to the correct orders of magnitude. This issue was revisited by Feldmann et 
al. who proposed to extract the mixing and decay constants on the quark flavor basis ■ This scheme can be easily 
extended to accommodate the mixing of ?7 C with 77 and 7/ from which the 772-77(77') mixing angles were extracted and 
turned out to be consistent with those from Refs. @, Q- 

Interestingly, the new data from BESIII for ip' — > 7P seem to suggest a deviation from the saturation assumption. 
It implies that some other mechanisms become important in ip' — > 7P, although they may not play a significant role 
in J/ip — > 7P. In this work, we shall show that the vector meson dominance (VMD) model is an ideal framework to 
make a coherent analysis of the 770-77(77') mixing effects and contributions from intermediate vector mesons. We shall 
show that the ip' — > jP is not saturated by the 77,3-77(77') mixing. Instead, one important mechanism that drives the 
difference between J/ip and ip' — > 7P and produces the observed patterns is the sizeable coupling of 0' — > J/ipP. 

As follows, we first give a brief introduction to the VMD model and lay out the correlated aspects of the 770-77(77') 
mixings in Sec. II. The detailed analysis, calculation results and discussions will then be presented in Sec. III. A brief 
summary will be given in Sec. IV. 
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II. VMD MODEL AND r/c-viv') MIXINGS 

In the VMD model (e.g. see review of Refs. [HI EH) the electromagnetic (EM) current can be decomposed into a 
sum of all neutral vector meson fields including both isospin-0 and isospin-1 components. The leading Vj* effective 
coupling can be written as: 

/V 7 = £^V^ (i) 

where V^(= p,w,(j>,J/ip f ...) denotes the vector meson field. The photon-vector-meson coupling constant eM v /fv 
can be extracted from the partial decay width IV_> e + e - . Neglecting the mass of electron and positron, we have 



e 

Jv~ 



3r, 



2a P 



(2) 



where p e is the electron three-vector momentum in the vector meson rest frame, and a e is the EM fine-structure 
constant. 

For the decays of J/0(0', -0(3770)) — > jP, the VMD contributing diagrams are illustrated in Fig. [1] This classi- 
fication is based on the photon producing mechanisms and related to the experimental measurements. For instance, 
Fig. [TJa) identifies such a process that the photon is connected to a hadronic vector meson fields. It requires a sum 
over all strong transitions of J/tp(ip', -0(3770)) — s- VP channels. 

The second process in Fig. [ljb) is via charmonium electromagnetic (EM) annihilations. Such a process generally 
has small contributions in comparison with the strong transitions. However, it is likely that the EM amplitudes 
may have significant effects in some exclusive decay channels. In recent series studies [12Ul5l] it shows that in the 
hadronic decays of J/ip(ip') — > VP, the short (via three gluon annihilation) and long-distance (Fig. [TJc)) transition 
amplitudes may have a destructive interfering mode that would efficiently reduce the strong transition amplitudes in 
some exclusive channels. As a consequence, the EM amplitudes may become compatible with the strong ones, and 
manifest themselves in experimental obscrvables. This issue is related to the so-called "pit puzzle" , which questions 
why the branching ratio fraction BR(tp' — > pit) j B R(J / 'tp — > pit) is so strongly suppressed in comparison with the 
pQCD expectation values E6H18l|. A review of this subject and some recent progresses on this problem can be found 
in the literature [H, [l§ HiP 

In the present work, our attention is to understand whether the data for J/0(0', -0(3770)) — > jP are consistent 
with those for J/ip(ip', 0(3770)) — > VP, and what drives the different radiative decay patterns between J/-0 and -0'. 
We shall adopt the available experimental measurements of J/-0(-0', 0(3770)) — > VP in the calculations of the VMD 
contributions. This means we need not worry about the detailed transition mechanisms for J/0(0' , 0(3770)) — > VP 
at this moment. Also, by adopting the experimental data for J/0(0', 0(3770)) — > VP, we need not consider the 
r/ — rf mixing processes since they have been contained in the data for J/0(0', 0(3770)) VP. 

It is worth noting in advance another feature with this classification of Fig. [1] Namely, transitions between vector 
charmonia may also contribute. For instance, 0' — > J/ipV will contribute to 0' — > 7?y. We will show later that this 
process is essential for understanding the radiative decay patterns for J/ip and 0' — > jP. 

Apart from the transitions via Fig. [TJ another important transition is via Fig. [2] which corresponds to the rj c -r)(r)') 
mixing due to the axial vector anomaly. Note that the process of Fig. [TJa) with an intermediate charmonium does 
not overlap with Fig. [5] at the hadronic level. In fact, it is interesting to note their correlated features: i) In both 
cases, the cc annihilate at short distances. In Fig. [Ha), the vector configuration of cc annihilates into a photon, 
i.e. cc in a relative S- wave with spin-1, while in Fig. [2] the pscudoscalar cc are in a relative S- wave but with spin-0, 
and then annihilates into gluons. ii) The process of Fig. [5] is through a typical magnetic dipole (Ml) transition 
of J/0(0', 0(3770)) — > jr] c , which can be regarded as a non-vector-resonance contribution in respect to the VMD 
scenario. 

With the Lagrangian of Eq. ([1]), the transition amplitude can be expressed as 



M 



VMD 



= fe ^-M^f^ ^vp) Hp^e^P^p^ , (3) 



V 7 



where g^vp denotes the coupling constants for the hadronic vertex of J/0(0', 0(3770)) — > VP, and will be determined 
by experimental data via 

M V p = g^vpHvD^uapP^pyey . (4) 
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FIG. 2: Schematic diagram for J/ip(ip' , ^(3770)) — > 777 and 777' via rj c -r](ri') mixing. 
We adopt an empirical form for the form factor [2ll - l23l |: 

where parameter /3 is in a range of 300 ~ 500 MeV. This form factor can be interpreted as the wavefunction overlap 
which would be suppressed in a large recoil momentum region for the final state particles [2ll - [23l |. The incovariant 
form factor can also be regarded reasonable in this case. The decay processes are treated in the cm. frame of 
the initial meson. Therefore, the anti-symmetric tensor structure of the interactions can always be reduced to a 
form of M^e^p ■ (jpy X €v), which explicitly depends on the three- vector momentum of the final state vector meson. 
Note that for the anti-symmetric tensor couplings all the contributions to the transition amplitude can be absorbed 
into the effective coupling form factor. Because of this, it is natural to expect that the form factor would contain 
information of meson wavefunction overlaps with an explicit three-vector-momentum dependence. In particular, a 
harmonic oscillator potential for the quark- antiquark system will lead to a form factor similar to Eq. ([5]). 

We shall determine the form factor parameter (3 combining the data for J/ijj(ip') — > VP and 7P. It will then be 
fixed and adopted for the calculations of other channels. In the transition of Fig. HJa), the vector meson will carry 
the momentum of the final state photon p 7 . 

The transition amplitudes of Fig. [2] can be expressed as 

» A mixinq \ t^u, v a B 

= Ap,/(p^^P;e> 7 a e? , (6) 

where \p Vc is the mixing angle between pseudoscalar P and r] c . It has been extracted in Ref. @, A I)I)c =— 4.6 x 10~ 3 
and Xn'n e = — 1-2 x 10~ 2 , which are also obtained by Ref. Q ■ It should be noted that in the above equation the coupling 
5i/>7r; c is extracted from the data for J/ip(ip') — > 777^,. The non-local effects from the off-shell ?7 C at the mass of 77(77') 
have been included in the mixing angles Q. In the second line, we define a reduced coupling Xpr k = Xp ric g^ nVc /P(p? f ) } 
which can be directly compared with the effective coupling (e/ fv)9ipvp in Eq. 

We do not include the 77^, mixings with the 77(77') in J/ip — > 777 and 777' since their mixing angles are relatively small. 
Nevertheless, the 77^, mixing effects will be further suppressed by the unknown but believe-to-be-small branching ratio 
for 77^ -> J/V>7- 
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Channels 


J/tp 


0' 


0(3770) 




- 


(1.30 ±0.10) x 10" 3 


< 2.8 x 10~ 4 


J/ipri 


- 


(3.28 ± 0.07)% 


(9 ±4) x 10 -4 


pit 


(1.69 ±0.15)% 


(3.2 ± 1.2) x 10 -5 


- 


p°7T° 


(5.6 ±0.7) x 10~ 3 




- 


UJ7V 


(4.5 ±0.5) x 10~ 4 


(2.1 ±0.6) x 10~ 5 


- 


0JT° 


< 6.4 x 10~ 6 


< 4 x 10~ 6 


- 


UJT) 


(1.74 ±0.20) x 10 -3 


< 1.1 x 10 -5 


- 


(fnj 


(7.5 ±0.8) x 10" 4 


(2.8tJ;°) x 10- 5 


(3.1 ±0.7) x 10~ 4 


p-q 


(1.93 ±0.23) x 10~ 4 


(2.2 ±0.6) x 10~ 5 




uirj 


(1.82 ±0.21) x 10 -4 


{3.2±H) x 10" 5 




<H 


(4.0 ±0.7) x 10" 4 


(3.1 ± 1.6) x 10" 5 




prl 


(1.05 ±0.18) x 10" 4 


(l.Qtll) x 10" 6 





TABLE I: Branching ratios for J/0(0', 0(3770)) -> VP from PDG 2010 Q. The dash "-" and dots ". . ." denote the forbidden 
and unavailable channels, respectively. 



III. NUMERICAL RESULTS 
A. Results from VMD 

In Table H the data for J/0, 0' and 0(3770) -> VP from PDG 2010 are listed. It shows that most of the 
light VP channels have been measured for J/0 and 0' hadronic decays. In contrast, most of the light VP channels 
for -0(3770) are below the experimental precision limit except for tprj. As mentioned earlier, the J/0(0') — > VP 
channels are correlated with the so-called "pw puzzle" in the literature. However, our attention in the present work is 
different. We shall use the experimental data for J/0(0', -0(3770)) — s- VP as an input to investigate the role played 
by the VMD mechanisms in the vector charmonium radiative decays. This treatment means that one need not be 
concerned about the detailed transition mechanisms for J/0(0', 0(3770)) — > VP at this moment since they all have 
been contained in the experimental data. We emphasize that this should not be a trivial starting point. Success of 
such a prescription would help us clarify two major processes in the charmonium radiative decays, i.e. the relative 
S'-wave cc annihilations would occur either via spin-1 or spin-0 configurations. 



In TablcQ] the branching ratios for 0' and 0(3770) —> J/ipV and J/07r° are also listed. As pointed out earlier, these 
channels are rather important for understanding the observed branching ratio patterns. The effective coupling g^pyp 
in Eqs. ^ and (jl]) is a scale-independent constant. The data in Table U will allow us to extract g^vp for different 
VP channels in association with the form factor parameter ft. The overall numerical study suggests that a smaller 
value of P = 0.3 GeV is favored. This is due to that in J/-0(0') — > jP, the intermediate vector mesons are in a highly 
virtual kinematic region. Part of the off-shell effects would be absorbed into the form factor parameter (3 as we adopt 
the V — > 7* couplings e//y which are determined by data for V — > e + e~ 0]. 

In Table |TT] we list the joint coupling constants (e/ fv)g^VP for different VP channels as a reference. These 
quantities are the corresponding scale-independent couplings in J/0(0') — > jP, and provide an immediate estimate 
of the relative strengths among those transitions amplitudes that involve different vector mesons. The form factor 
J^p^) = e~ p i/ 8 ^ with (3 = 0.3 GeV will lead to an overall suppression to the vertices. In the light VP sector, the 
strong /9°7 coupling accounts for the relatively large contributions from the p° mediated transitions. 

In the vector-charmonium-mediated channels, the non-negligible coupling of g^,'j/^ v implies a non-vanishing cou- 
pling of g^,'j/^ n ', although the decay of 0' — > J/iprj' is prohibited by the phase space. The influence of 0' — > J/0r/ in 
0' — > 77/ should not be neglected and must be included in the amplitude. As we know, the rj and 7/ can be expressed 
as mixtures of quark flavor singlets: 

i] = cos apnn — sin apss , 

r/' = sin apnn + cos apss , (7) 
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VP 


(e/ 'fv)g.j/i,vp 


(e/fv)g-4,>vp 


( e / fv)g^(3770)VP 


J/lP-R° 


- 


4.02 x 10~ 4 


< 9.73 x 10" 4 




- 


6.25 x 10" 3 


2.74 x 10" 3 


J/W 


- 


3.01 x 10 -2 


1.32 x 10 -2 




2.40 x 10~ 4 


- 


- 


tp'r] 


3.74 x 10" 3 


- 


- 




1.80 x 10 -2 


- 


- 


V>(3770)tt° 


< 3.22 x 10~ 4 


- 


- 


■0(3770)77 


9.09 x 10" 4 


- 


- 


^(3770)?/ 


4.37 x 10~ 3 


- 


- 


p%° 


2.83 x 10" 3 


6.69 x 10" 4 






2.35 x 10~ 4 


2.73 x 10" 4 






< 2.51 x 10 -5 


< 1.02 x 10~ 4 






3.97 x 10" 4 


< 1.67 x 10" 4 






2.72 x 10~ 4 


2.69 x 10" 4 


1.02 x 10 -2 


P7? 


4.52 x 10 -4 


8.10 x 10~ 4 




uirj' 


9.54 x 10" 5 


2.02 x 10" 4 






1.48 x 10~ 4 


1.20 x 10 -4 




pv' 


2.48 x 10~ 4 


5.34 x 10~ 4 





TABLE II: Effective couplings j^g^vp (in unit of GeV 1 ) for J/t/j(t//, 0(3770)) — > 7P extracted from the intermediate VP 
channels. Note that the form factor ^(p 2 ) is not included. The dash "-" and dots ". . ." denote the forbidden and unavailable 
channels, respectively. 



where ap = arctan-s/2 + Op, and Op ~ —11.7° is the SU(3) flavor singlet and octet mixing angle. Thus, we have 

(\f2smap + Rcosap \ 
u~- ' ( ' 
V2cosap — Rsmap J 

where R describes the SU(3) flavor symmetry breaking. In general, R = f-^/ fi< — 0.838 is commonly adopted for the 
relative production strength of an ss to qq. The above relation is based on the qq and ss mixing scheme [HI, I24rl27j 
and does not include a possible glueball component. If one extends the 77-77' mixing to accommodate the glueball Q , 
the coupling of g^'j/^ 1 can be expressed as 

( V2X n , + RY V > + GZ„, \ , n , 

where parameter G denotes the relative strength of producing the pseudoscalar glueball Q to a light qq component. 
The general flavor wavcfunctions for 77 and 77' are 

77 = Xr,nn + Y tl ss + Z^Q , 

77' = X,fnn + Y,fss + Z,yG , (10) 



for which different model solutions can be found in the literature [12l. l24l - l26l .l28l]. Generally speaking, the introduction 
of the glueball component will introduce new parameters. Taking into account that the glueball components within 
the 77 and 77' are rather small, and Eq. ([7]) is well established to leading accuracy, we neglect the possible glueball 
mixing effects in the present analysis. 

We adopt the same on-shell couplings of gj/^'p as those extracted in ip' — > J/tpP since the kinematics for these 
two processes are similar to each other. Namely, we neglect the off-shell effects with the couplings of gj/^>p in 
contrast with g^/j/^p- 

As listed in Table [TTl it shows that the charmonium poles are one of the most important contributing sources to 
the J/tfj(ip') — > 777 and 777', which seems to be slightly out of expectation and has not been addressed before. This 
feature is explicit for the ip' decays since the decay of tp' — > J/ipV is experimentally accessible. In contrast, other 
VP channels' contributions to 7P are rather small due to their relatively small branching ratios. Similar phenomena 
appear in J/tp — > VP except that the sizeable branching ratio for J/ip — > pn would also make the pix channel an 
important contributor to the 7P amplitude. 
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WV'^7P) 


7?? 
777' 


2.10 x icn 2 

3.66 x 10" 2 


5.12 x 10~ 3 
8.91 x 10~ 3 



TABLE III: Reduced effective couplings (in unit of GeV 1 ) from the r] c -ri(ri') mixings. 





J/V -»• 7P 


ip' -» 7P 


7 P 


VMD 


r/ c mixing 


VMD 


T) c mixing 


77T° 
777 
777' 


(1.64 ~ 2.04) x 10" 5 
(0.060 ~ 0.063) x 10~ 3 
(1.04 ~ 1.05) x 10" 3 


0.61 x 10~ 3 
3.50 x 10" 3 


(0.66 ~ 1.15) x 10~ 7 
(3.33 ~ 3.61) x 10~ 6 
(0.58 ~ 0.61) x 10 -4 


1.62 x 10~ 6 
0.096 x 10" 4 



TABLE IV: Branching ratios for J/ip(ip') — > 777 and 777' given by the VMD mechanisms and 770-77(77') mixings, respectively. 



B. Results from 770-77(7/) mixings 

In Tabic ILnl we list the effective couplings derived from the 77,^-77(77') mixings 0|. These values can be directly 
compared with (e/ fv)9i>VP listed in Table ILT1 It shows that in J/ip — > 777 and 777', the axial-anomaly-driving mixing 
contributions turn out to be more predominant than the VMD, while in ip' — > 7P the most important contribution is 
from the J/ip pole. 

We list the individual branching ratios given by the VMD and r) c -T)(j)') mixings in Table UVl as a comparison. Indeed, 
it shows that the mixing contributions have nearly saturated the branching ratios in J/ip — > 777 and 777'. However, 
the situation changes in ip' decays where the VMD mechanisms become more important. An interesting feature is 
that one in principle needs both to give an overall account of the measured branching ratios. 

Note that in Table IIV1 the ranges of uncertainties for the VMD results are given by the experimental errors in 
Table H 



C. Discussions 



To compare with the experimental measurements, we need to add the VMD and 77,3-77(77') mixing amplitudes to 
each other coherently. Taking the advantage of the unique Lorentz structure of the VVP coupling, we can express 
the total transition amplitude as follows, 

M tot = M™ D + e i5 M™ 3 , (11) 

where 5 is introduced to take into account possible phase differences between these two amplitudes. In the transition 
processes that we are interested in here, such a phase ambiguity seems inevitable due to a important role played by 
hadronic transition mechanisms. Since several different hadronic level amplitudes are involved in Ai^p ID , it is not 

a necessity that Ai^p ID and _/vf™ xm 9 should share the same phase angle for different pseudoscalar channels. We 
expect that the experimental data 0, 0] would provide a constraint on it. 

In Fig. [3J we plot the (^-dependence of the branching ratios in comparison with the PDG2010 averages Q and new 
experimental data from BESIII 0|. It shows that in the two decays, J/ip — > jrj and ip' — > 777', the transition amplitudes 
of the VMD and 77^77(77') mixings are well in phase. In contrast, they seem to be out of phase in ip' — > 777, although 
the experimental uncertainties are quite large. The central value of the data can be best accounted for at S ~ 140° 
or 220°. More complex phases appear in J /ip — > 777', although the dominant contributions are from the axial gluonic 
anomaly. In this case, the phase angle S = 80° or 280° are favored. It should be mentioned that in a recent paper by 
BESIII Hi|, a smaller branching ratio for J/ip -> 777' is reported, i.e. BR(J/ip -> 777') = (4.86 ± 0.03 ± 0.24) x 10~ 3 . 
This value is consistent with the PDG2010 average, and would favor S ~ 90° or 270°. 

In Table fVl we list the coherent results for the branching ratios BR(J/ip —> 7P) and BR(ip' — > 7P) in comparison 
with the data again 0,01. The phase angles are fixed as shown in Fig. [3] with the best description of the central 
value of the data. Again, the theoretical uncertainties due to adopting the data for J/ip(ip') —> VP are included. 




FIG. 3: The phase angle (5-dependence of branching ratios for J/ip — s- 777''' ((a) and (b)) and ip' — ► 777''' ((c) and (d)). The 
experimental data with uncertainties are shown as the straight bands, while the theoretical results are shown by the curvilinear 
bands. The theory uncertainties are given by the uncertainties of the data for J/0(0') — > VP as listed in Table|T] 





j/tp -> yP 


0' -> 7P 


0(3770) ->■ 7P 


IP 


Experiment 


Theory 


Experiment 


Theory 


Experiment 


Theory 


77T° 

IV 

W 


(3.49±g;») x 10- 5 
(1.104 ±0.034) x 10~ 3 
(5.28 ±0.15) x 10" 3 


(1.64 ~ 2.04) x 10" 5 
(1.05 ~ 1.06) x 10" 3 
5.20 ~ 5.22 x 10~ 3 


(1.58 ±0.40 ±0.13) x 10" 6 
(1.38 ±0.48 ±0.09) x 10" 6 
(1.26 ±0.03 ±0.08) x 10" 4 


(0.66 ~ 1.15) x 10 -7 
(1.39 ~ 1.53) x 10~ 6 
(1.14 ~ 1.19) x 10" 4 


< 2 x 10" 4 

< 1.5 x 10~ 4 

< 1.8 x 10~ 4 


3.25 x 10" 9 
7.95 x 10" 7 
1.64 x 10" 6 



TABLE V: Calculated branching ratios for J/0(0', 0(3770)) -> ■yP based on the VMD model. Experimental data from PDG 
2010 H for J/ip and V(3770) decays and from BESIII [3 for ip 1 decays are included as a comparison. The phase angles are 
fixed in such a way that the theoretical results can best describe the central values of the experimental data. 



We also include the ?^(3770) — > 7P as a prediction of the VMD mechanism. The predicted branching ratios are all 
small. The rj c mixing contributions are not included here due to lack of data. Also, most of the light vector meson 
contributions to the ^(3770) radiative decays are rather small and unavailable. Thus, the predicted branching ratios 
are actually given by the J/ip pole in the VMD model. Given the same statistics for ■0(3770) as the ijj' from BESIII, 
the accessible channel would be -0(3770) — > 777'. Experimental examination of the predicted pattern in Table fVl would 
be an interesting test of the VMD mechanisms proposed in this work. 

In general, the results fit the observed branching ratio pattern very well, except that the branching ratio for 
0' — > 77T° seems to have some discrepancies. It might be a sign that other non-VMD mechanisms may also play a 
role. For J/ip —> 771" , the dominance of J/0 — > p°Tr° can naturally account for the data. It should be mentioned that 
Ref. [13] also confirms the VMD contributions via the p°ir° channel to J/0> — > jit . 
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Our investigation suggests the importance of a coherent treatment for the VMD mechanism and rj c -ri(r]') mixings. 
Note that the charmonium pole contribution has not been included by the previous studies 0-0; EH • Meanwhile, 
an understanding of why the VMD and axial gluonic anomaly mechanisms play different roles in J/ip and ip' decays 
would be essentially important. The following points may help to clarify this question: 

i) As mentioned earlier, there are some interesting correspondences between the axial gluonic anomaly and VMD 
in this case. In the axial gluonic anomaly transitions the cc annihilate into gluon fields at short distances in a relative 
S'-wave and spin-0, which induces mixings with the Goldstone boson rj and SU(3) flavor singlet 77'. The photon 
radiation can be regarded as from non-vector-resonance Ml transitions. In the VMD transitions via the charmonium 
state, the cc also annihilate at short distances in a relative S'-wave, but with spin-1. In this case, the annihilated cc 
couple to a photon, and radiate two soft gluons which can couple to pseudoscalar states. 

ii) The difference between those two mechanisms can be well-understood quantum mechanically. For J/ip — > 777^, 
the VMD transitions via ip' pole is relative suppressed by the tp'^y coupling since as the first radial excited state the 
wavefunction of ip 1 at the origin is smaller than that of J /ip. In contrast, the axial-gluonic- anomaly-driving T] c -r](r)') 
mixings will occur via J/ip — > jr] c — > 77/'', where the first step is a typical EM Ml transition between two ground 
charmonium states. It is allowed by the quantum transition selection rule at leading order. 

The situation changes in ip' — > 777W. On the one hand, the VMD transition will be dominated by the J/ip pole, 
which will be coupled to the EM field. On the other hand, the axial gluonic anomaly transitions via the ip'-rj c Ml 
transition will be suppressed by the quantum transition selection rule at leading order. For the 77^,-mcdiatcd transition, 
the r]' c mixings with the r\ and 77' will then be suppressed . 

The above qualitative argument explains why the VMD mechanism and axial gluonic anomaly play different roles 
in J/ip and ip' decays, respectively, as manifested by the calculation. In particular, it shows that both mechanisms 
are crucial for our understanding of the observed branching ratio patterns. 

The successful account of the observed branching ratio patterns for J/ip(ip') — > 7P in the VMD model has an 
important implication of the hadronic decay mechanisms for J/ip(ip') — > VP. It shows that the "puzzling" radiative 
decay patterns in J /ip{ip') — ► 7P have direct connections with the hadronic decay mechanisms, i.e. J/ip(ip') —> VP, 
instead of some other abnormal processes. As a consequence, it will guide our further investigations of the transitions 
of J /ip(ip') — > VP, and impose constraints on processes such as illustrated by Fig. [1] For instance, the hadronic 
part of Fig. [He) is found to be an important non-perturbative QCD mechanism that contributes predominantly in 
ip' — > J/ipTj and J/ipir° [3l|, [S^]. As pointed out recently in a series of papers on the subject of non-perturbative 
transition mechanisms in charmonium decays (hi . I20I I3ll - l35l ] , such intermediate meson loop transitions would be an 
natural mechanism for evading the pQCD hclicity selection rule and explaining the "pit puzzle" in J/ip(ip') — > VP. 

IV. SUMMARY 

In brief, with the available data for J/ip(ip') — > VP, we show that the VMD model is still useful for our understanding 
of the newly measured branching ratios for J/ip(ip') — > 7P in association with the r/c-viv') mixings via the axial 
gluonic anomaly. Importance of such a contribution has not been recognized before. In particular, we stress that 
the intermediate vector charmonia can have significant contributions via e.g. ip' — >• J/iprj — > 777. We show that these 
two mechanisms behave differently in J/ip and ip' — > 7P, and can be understood by state transition selection rules. 
We also emphasize that the consistency between J/ip (ip 1 ) — > 7P and VP demonstrated in this work would impose 
important constraints on the non-pQCD mechanisms in J /ip{ip') VP. It would be useful for our final understanding 
of the long-standing "/J7r puzzle" in J/ip(ip') — > VP. 
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